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DIFFUSION  OF  THE  SELF  MAGNETIC  FIELDS  OF  AN  ELECTRON  BEAM 
THROUGH  A  RESISTIVE  TOROIDAL  CHAMBER 

I.  Introduction 

Extensive  work  has  been  done  in  the  past  on  the  diffusion  of  external 
■agnetic  fields  into  hollov  circular  cylindrical  conductors  of  infinite 
length^  ^  and  also  on  the  diffusion  of  the  self  magnetic  field  of  beams  out 
of  hollov  cylinders.  ~  The  beam  self  magnetic  field  diffusion  studies 
have  also  considered  the  effect  of  the  diffusing  fields  on  the  beam 

6—8 

dynamics  and  have  furnished  Interesting  results  on  the  beam  stability 

9 

and  beam  trapping.  However,  these  studies  are  linear  and  the  expressions 
for  the  fields  are  valid  only  near  the  axis  of  the  cylinder. 

Following  the  installation  of  strong  focusing  windings  in  the  NRL 
modified  betatron  accelerator  it  is  routinely  observed  that  for  several 
combinations  of  injection  parameters  the  injected  beam  consistently  spirals 
from  the  injection  position  to  the  magnetic  minor  axis  and  is  trapped. 
Attempts  to  explain  this  interesting  phenomenon  using  the  existing  linear 
resistive  model  have  been  unsuccessful.^^  The  decay  rate  predicted  by  the 

9 

linear  theory  for  the  parameters  of  the  experiment  is  at  least  10-20  times 

longer  than  that  observed  in  the  experiment,  even  when  wake  field  effects 

12 

are  taken  into  account. 

9  12 

In  contrast  to  the  analysis  '  that  assumes  the  bean  to  be 
near  the  minor  axis,  the  beam  in  the  experiment  during  injection  is  at 
least  temporarily  near  the  wall.  In  addition,  the  geometry  of  the  NRL 
device  is  toroidal^^  and  not  cylindrical  and,  therefore,  there  are 
additional  characteristic  times^  that  may  modify  the  diffusion  process. 

The  present  work  extends  the  results  of  the  linear  theory.  The 
expressions  for  the  diffusing  fields  are  valid  not  only  near  the  axis  but 
almost  over  the  entire  cross-section  of  the  chamber  and  toroidal  effects 
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are- included  to  lowest  order.  The  results  have  been  obtained  under  the 
following  assumptions.  First,  it  is  assumed  that  the  toraidal  vacuum 
chamber  has  a  small  aspect  ratio,  i.e.,  the  results  are  valid  provided  that 
the  ratio  of  the  minor  to  the  major  radius  is  much  less  than  unity  and  the 
radial  distances  of  the  observation  point  from  the  minor  axis  is 
considerably  smaller  than  the  major  radius  of  the  torus.  Second,  since  the 
results  are  confined  in  the  vicinity  of  the  toroidal  chamber,  propagation 
effects  do  not  play  any  role,  i.e.,  the  displacement  current  is  omitted  in 
Maxwell's  equal'ions.  Third,  in  order  to  obtain  tracktable  analytical 
results,  it  is  assumed  that  the  conducting  wall  is  thin,  i.e.,  its 
thickness  is  much  smaller  than  the  minor  radius  of  the  torus.  In  this 
case,  the  analytical  results  are  not  valid  very  near  the  conducting  wall. 

In  the  limit  when  the  ratio  of  the  wall  thickness  to  the  minor  radius  of 
the  torus  tends  to  zero,  i.e.,  for  a  toroidal  conducting  shell,  the 
analytical  results  are  exact  everywhere  inside  of  the  toroidal  vessel. 
Finally,  the  analytical  results  on  the  beam  dynamics  are  further  simplified 
under  the  assumption  that  the  current  ring  moves  slowly  in  comparison  to 
the  fastest  of  the  characteristic  times  that  dictate  the  diffusion  process. 

Under  the  assumptions  mentioned  above,  it  is  found  that  there  are 
three  characteristic  times  with  which  the  magnetic  field  leaks  out  of  a 

resistive  torus,  when  a  current  ring  turns  on  at  t  -  0  inside  the  torus. 

2  2 

The  shortest  is  the  "plane"  or  "fast"  diffusion  time  Tpp  -  u^aCb-a)  /n  , 
where  a  is  the  wall  conductivity,  a  and  b  are  the  inner  and  outer  minor 
radii  of  the  conducting  wall  and  is  the  permeability  of  the  vacuum.  The 
terms  associated  with  are  responsible  for  the  electric  field  to  be  zero 
at  tsO  outside  the  torus  since  no  leakage  has  occurred  as  yet.  The 
"cylinder"  diffusion  time  tp  «  u^a(b-a)a/2  together  with  Xpp  determines  the 
speed  with  which  the  fields  penetrate  the  wall  chamber  so  that  the  images 
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of  these  fields  gradually  disappear.  Finally,  the  "loop"  diffusion  tine 
too  *  ~  vhere  r^  is  the  najor  radius  of  the  torus, 

determines  the  speed  with  which  the  fields  diffuse  into  the  hole  of  the 
doughnut  and  is  responsible  for  the  gradual  disapi^arance  of  the  wall 
current.  It  turns  out  that  the  radial  component  of  the  self  magnetic 
field,  which  is  responsible  for  the  beam  trapping  in  the  MBA,  is 
independent  of  and,  therefore,  the  "loop"  diffusion  tine  does  not  play 
any  role  in  the  resistive  trapping  of  the  beam.  The  expressions  of  the 
fields  predicted  by  the  present  work  have  been  used  to  compute  the  beam 
centroid  orbit  and  several  other  trapping  parameters  measured  in  recent 

iq 

detailed  beam  trapping  studies  in  the  NRL  modified  betatron  accelerator. 
The  shape  of  the  computed  orbits  is  very  similar  to  those  observed  in  the 
experiment . 

In  Section  II,  the  diffusion  problem  for  a  current  ring  inside  a 
toroidal  conductor  is  formulated.  In  Section  III,  the  vector  potential  for 
a  current  ring  in  the  absence  of  the  conductor  is  derived.  This  is  the 
particular  solution  of  the  problem.  In  Section  IV,  the  homogeneous 
solution  inside  and  outside  the  torus  (but  not  inside  the  conductor)  is 
derived  in  toroidal  geometry.  The  toroidal  geometry  removes  the  ambi{(uity 
on  the  value  of  certain  constants  associated  with  the  logarithmic 
dependence  of  the  solutions.  In  Section  V,  the  initial  conditions  are 
established  that  determine  the  time-dependent  arbitrary  parameters  in  the 
homogeneous  solutions.  In  Section  VI,  the  solution  inside  the  conductor  is 
derived  and  the  boundary  conditions  are  applied.  The  vector  potential  is 
computed  in  the  three  regions  Inside  and  outside  the  torus  and  inside  the 
conductor.  In  Section  VII,  approximate  results  are  obtained  under  the 
assumption  of  a  thin  conducting  wall.  Section  VIII  contains  exact  analytic 
results  for  the  shell  model  and  Section  IX  provides  a  summary  of  the  main 
results  and  lists  the  most  important  conclusions  of  the  present  study. 
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II.  -PorBulation  of  the  Diffusion  Problee 

The  configuration  and  system  of  coordinates  is  shovn  in  Pig.  1.  The 
toroidal  chamber  has  a  major  radius  r^,  an  inner  and  outer  minor  radius  a 
and  b,  respectively,  and  conductivity  a.  In  the  presence  of  an  external 
driver,  namely,  a  current  ring  which  is  axisymmetric,  time-dependent  and  is 
located  inside  the  torus,  the  magnetic  vector  potential  in  that  region  is 
determined  by  the  equation 

^  X  ^  X  -  uj,  (1) 

where  the  vector  potential  t  has  only  one  nonzero  component  that 

depends  only  on  the  cylindrical  components  (r,z)  and  on  tine.  The  current 

density  1  of  the  current  ring  has  only  one  nonzero  component  which  is 

2 

equal  to  ■  I^/nr^  inside  the  ring  and  zero  outside  it.  Here,  I^  is  the 
o  c  c  ^  c 

ring  current  and  r^  is  its  minor  radius. 

The  magnetic  field  inside  the  conductor  is  determined  from  the  vector 
potential  that  is  described  by  the  diffusion  equation 

^  X  ^  X  -  -  Uq  "  It^' 

where,  again,  has  only  one  component  which  depends  on  (r,z)  and 

on  time.  Finally,  the  vector  potential  outside  the  torus  is  determined  by 
the  homogeneous  equation 

^  X  ^  X  -  0  (3) 

where  the  component  A^^*  of  depends  on  (r,z)  and  on  time.  The 
magnetic  and  electric  field  components,  in  each  region,  are  given  by 
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where  is  one  of  the  components  defined  above,  depending  on  the  region  of 
interest.  Notice  that  Eq.  (2)  is  identical  to  Ampere's  law  combined  with 
Ohm's  law  inside  the  conductor. 

It  is  convenient  to  express  the  vector  potential  and  the  fields  in 
terms  of  the  local  cylindrical  coordinates  (p, ♦)  which  are  related  to  the 
global  cylindrical  coordinates  (r,z)  by: 


r  -  r^  +  pcos^, 


z  >  psin#. 

Then,  in  the  region  inside  the  conductor,  Eq.  (2)  reduces  to 


1  a_ 

p  ap  **  ap 


.2. con 
1  ® 

+  T  “2 — 


aA|°"  .  aA|°" 

JT"  -  p  a^^ 

r^  +  pcos4 


(5a) 

(5b) 


A  con 

2 

(r^  +  pcos4) 


-  + 


at  » 


(6) 


after  taking  account  of  the  fact  that  Aq°"  is  independent  of  the  toroidal 
angle  9.  Also,  in  the  local  coordinate  system  Eqs.  (4a),  (4b)  are  replaced 
by 


B 

p 


1  sin4  . 

"  p  a4  r^  +  pcos4  ^9’ 


(7a) 
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(7b) 


R  ^  cos*  . 

*♦  "  3p  *  r^+pcos4  ^0’ 


while  Eq.  (4c)  renains  the  saee. 


At  the  surface  of  the  toroidal  conductor»  i.e.,  at  p  ■  a  and  p  -  b, 
the  tangential  components  of  the  electric  and  magnetic  fields  are 
continuous.  Therefore,  in  the  local  coordinate  system,  the  boundary 
conditions  are: 


Eq"  (p  -  a,  ♦,t)  -  (p  -  a,  *,1), 


C'  (P  -  b,  ♦,t)  -  eJ°"  (p  «  b,  ♦,!), 


bJ"  (p  -  a,  ♦,t)  .  bJ°"  (p  -  a,  ♦,!), 


bJ"*  (p  -  b,  i,t)  -  bJ°"  (p  .  b,  ♦,!). 


Since  the  vector  potential  is  sero  at  t  >  0,  the  first  two  boundary 
conditions  can  also  be  expressed  as 


aJ"  (p  -  a,  ♦,!)  -  A®""  (p  .  a,  ♦,t), 


aJ"'  (p  -  b,  ♦,t)  -  A^°"  (p  .  b,  ♦,!). 


Thus,  the  diffusion  fields  in  the  three  regions  inside  and  outside  the 
torus  and  inside  the  conductor  are  determined  by  the  solutions  of  Eqs.  (1), 
(3)  and  (6)  with  the  boundary  conditions  given  by  Eqs.  (9a),  (9b),  (8c)  and 
(8d)  on  the  inner  and  outer  surface  of  the  toroidal  conductor. 
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Ill i  Vector  Potential  of  a  Current  Ring 

In  this  section,  an  approximate  expression  for  the  vector  potential  of 
a  current  ring  external  driver  is  obtained.  This  is  the  particular 
solution  of  Eq.  (1).  For  that  purpose,  the  toroidal  coordinate  systee 
associated  with  the  current  ring  is  used.  Toroidal  coordinates  are  eost 
appropriate  for  a  toroidal  conductor  and  their  significance  vill  becoee 
evident  in  the  next  section  vhen  the  arbitrariness  as  to  the  value  of 
certain  constants  in  the  solution  is  reeoved. 

The  global  cylindrical  coordinates  (r,z)  are  related  to  the  ring 
toroidal  coordinates  by: 


b  sinhn* 

“c  coshh'  -  cos  V ' 


(10a) 


I  sinC* 

'c  coshh'  -  cos  * 


(10b) 


where  b^  is  a  constant.  These  relations  can  be  easily  inverted,  naeelyt 


(r  -  b^)^  4  2^ 
(r  +  b^)2  +  2^’ 


(11a) 


r^'cos  V  -  j  (i 

(lib) 

sin  5'  -  Ip 

(X  . 

(11c) 

According  to  Eq.  (11a),  vhen  h'  is  fixed,  the  coordinates  (r,z)  describe  a 

circle  whose  radius  is  b^/sinhq' .  If  for  q'  ■  this  circle  coincides 

with  the  current  ring  surface  whose  major  radius  is  R^,  then  it  is 

2  2  1/2 

straightforward  to  show  that  b^  -  [R^  -  r^]  .  The  points  (h',^')  outside 
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the  ring  are  detenined  by  the  inequality  h'  <  idiile  the  points  inside 
the  ring  are  deterained  by  h'  > 

In  the  absence  of  the  toroidal  conductor,  the  solution  of  Eq.  (1)  for 

2 

a  current  ring  with  constant  current  density  ■  l^/nx^  Is  equal  to 


.ext  *^0  c  cos(6'  -  8") 

®  All^r  2  .  in  I 

®  V 


where  x'  and  x”  are  the  observation  point  and  a  point  inside  the  ring, 
respectively,,  and  the  integration  is  over  the  volune  V  of  the  ring.  The 
Green's  function  of  |x'  -  x"|~^,  in  toroidal  geoaetry  is  equal  to 

— T — ~ (cosh  h'  -  cos  5')^^2(cosh  n"  -  cos 
|x'  -  I  "‘‘c 


^  n  ^(a-n+j) 

)  - f-  cosn(e'-e»)cosa(t'-t") 

6.0  r(-n4) 


Ipn 


j/2  (cosh  n')  qJJ  _ 


'  (cosh  q")  0  (cosh  r\ 

a  -  1/2  a  -  1/2 


-);  r  >  h'  \ 
');  h"  <  n' 


where  e  «  1,  e  -  2  when  a  -  1,  2,  3,  ...,  r(x)  is  the  gaaaa  function,  and 
o  a 

P"  , (cosh  n)»  Q"  1/0  (cosh  n)  are  the  associated  Legendre  functions 
n  -  i/z  B  -  1/a 

of  the  first  and  second  kind,  respectively.  Without  giving  the  details  of 

the  calculation,  substitution  of  Eq.  (13)  into  Eq.  (12)  leads  to 

the  following  expression  of  the  vector  potential  of  the  current  ring: 


-  b^(cosh  n'  -  cos  ^ 


oi  1/9  (cosh  h')  cos  a  V 

B  B  B  l/A 
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b^(cosh  ri'  -  cos  y -  1/2  * 


(Ua) 


Inside  the  current  ring,  i.e.,  when  h'  >  h'  f  and 

c 


^0*^  -  b^(cosh  h'  -  cos  ®b**b*^^b  -  1/2  ®°®  ■  » 


(14b) 


outside  the  current  ring,  i.e.,  when  h'  <  h'^..  where 


4/2y  I  b  1  r’’'  1  9  .1  I. 

3„2^  2  _  1  J  **■  -  1/2  -  1^2  sinh  h"’ 


(15a) 


A/Iy  I^b 


I  J  Jj 

rP  ^  «i  - 1/2  "">0^  - 1/2  <“»>'  ""> 

■  -  r  J 


(15b) 


4/?y  I_b 


I  b  r 

rP  y-T  <>i  - 1/2  <“*•■  ’’">“2  - 1/2  <■=«•' 

®  -  r  J 


dn" 

sinh  H"’ 


(15c) 

and  h'  -  tn  [(R  +  r  +  b  )/(R  +  r  -  b  )].  In  the  derivation  of  Eqs. 

w  c  c  c  c  c  c 

(14a),  (14b)  from  Eq.  (12),  use  was  made  of  the  identity 


- SOSSLi: - -  dr  -  o-  Q„^  1/2  <*'°®** 

(cosh  n"  -  cosr)  3  sinh^^n"  ~ 
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The  toroidal  functions  P”  _  (cosh  h)  «nd  Q®  _  *1) 

appearing  in  Bqs.  (13),  (14a),  (14b),  (15a),  (15b)  and  (15c)  are  given  by 
the  following  exact  expressions: 


Pj  .  (cosh  n)  -  1,  ‘“‘"h  «”• 

It  Hn  -  n  +  j*/ 


(m  .  n  -  7)n 


m-l  (n  +  ^)  (n  -  a  +  x) 


*  (1  - 


««o> 


I 

s-0 


s!  (I  -  a). 


_s  2sll 


(_l)V^^r(m+n4)  „ 


-  (B+n+j)n 


•  (n+j)  (m+n+j) 


Z  — ^  [‘"  <*•’’>  *  "s  *  V. 


V  -  V 

n+s  B+n+s 


]. 


SaO 


(17a) 


(-l)"2“n^^^r(in+n+|>) 


-  1/2  - nm*T) 


(sinh  h) 


(m-«-n+j)r) 
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I 

SmO 


(n+i)  (B+n+y) 

*^5  s  -  2sh 

Sl(B+l)g  ® 


where 

(a)  ■  a(a+l)(a+2) . • • (a+s— 1),  (n)Q  *  1» 

S  " 


(17b) 


(18a) 


(18b) 


'n  -  ^  2  1' 


-  0. 

o 


(18c) 


For  B-O,  the  first  term  in  Eq.  (17a)  is  omitted,  since  •  1  for  b  >  n, 
and  -  0  for  man.  The  expressions  above  are  appropriate  for  the 
region  inside  as  veil  as  outside  the  torus  but  on  its  vicinity. 

Up  to  this  point,  the  results  given  above  are  exact.  In  the 
following,  an  approximate  expression  of  vill  be  obtained  froB  Eqs. 
(14a)  and  (14b)  under  the  assumption  of  a  small  aspect  ratio  ot  the 

“^c 

current  ring  by  keeping  terms  up  to  order  e  .  Notice  that,  vhen  r^/R^ 
«  1,  then  to  first  order  in  the  aspect  ratio,  or  to  first  order  in  the 


toroidal  corrections,  ve  have  the  approximate  relations:  b^  e  R^,  e  ^ 

e  rg/2R^  and  e”*^  **  p'/2R^.  Here,  (p',  ♦')  are  the  local  cylindrical 
coordinates  vith  respect  to  the  ring  position,  i.e.. 


r  ■  R^  +  p'  cos  ♦' , 


(19a) 


z  -  p'  sin  4'. 


(19b) 


Making  use  of  the  identity 


(cosh  h'  -  cos  5') 


,,1/2 


Z 


(n')  COSB  t', 


where 


D^(h') 


[o_  (cosh  V)  -  prjlT  *'’>]’ 


(21a) 


11 


■  "  vi/2<“»'' "'))],  <2“) 

and  also  of  Eq.  (17b),  a  straightforward  calculation  leads  to  the 
approximate  expression,  to  first  order  in  toroidal  corrections, 

m 

b^(cosh  t\'  -  cos?')^^^  ^  ^n  -  1/2 

rnmo 

“-^*>0  [ar  -  (ar  -  3  a  -S  e-»^'cos  (22) 

The  number  of  terms  kept  on  the  right  hand  side  of  Bq.  (22)  was  determined 

by  the  fact  that  the  quantities  b  af^^  are  of  zero  order  in  toroidal 

corrections,  as  is  indicated  by  Eqs.  (15a),  (17a)  and  (17b).  Similarly, 

the  quantities  b^  b***  and  b^  c^^  are  of  order  (b^/r^)^  g-2(a+l)ir|' 

C  HI  c  ■  c  c 

2  -2(m+l)n' 

(b^/r^)  e  and,  consequently,  we  have  the  approximate  expression 

m 

b  (cosh  h'  -  cosS')^^^  }  , ,,  (cosh  h')  cosm  C' 

c  mm  m  -  i// 

m>o 

■  -  H  ‘'c  [C'  <*»"'>  -  2) 

-  b®**  (en  (46*^')  -  2)  e"^'  cos  I'  -  2bJ*^  e^'  cos  K’  ],  (23) 

to  first  order  in  toroidal  corrections.  The  same  relation  above  holds  for 

the  quantities  b  c_  . 

cm 

From  Eqs.  (14a),  (14b),  (19a),  (19b),  (22)  and  (23),  it  follows  that 
in  the  local  coordinate  system  (p',  f')  of  the  current  ring  and  to  first 
order  in  ve  have  inside  the  ring,  i.e.,  when  p'  <  r^: 


—  K 
2/1  ^ 


ext 

1  1 

f  ext  -  ^ext)  d'cos4' 

*0 

-  ?  1 

1*0  -  5  *1  J  rT^J 

r  8R  a 

f 

2r  \ 

®o  1 

tn  -p-  -  2J  + 

V 

_o _ J 

»'cos4' 

■ren!!^.2l 

1 

-  4cj 

2 

r"  p'  J 

,  r^^cos+'  ' 

d^cosA* _ c 


R. 


ReP' 


(24a) 


while  outside  the  ring,  i.e.,  when  p'  ^  r^,  we  have: 


•  ext 

ftn  2I 

C'  r 
*  ° 

a 

2n\ 

^o 

r"  p'  ■ 

--r  [ 

cos»-1 


>r' 

2 


“r'  .'COS*'  ,.ext  V“*' 


llfMiL  (t„  -  2)  -  4bf 


(24b) 


Next,  we  need  the  approximate  expressions,  to  first  order  in 
for  the  quantities  *c**n*^'  m  ■  0,  1.  Prom  Eqs. 

(15a),  (15b),  (15c),  (17a),  and  (17b),  it  is  easy  to  show  that 


R  ^ext 

R.a  a  —  n 
CO 


“c  3  f,  “'•c  3ll 


(25a) 


Iw- 


(25b) 


if]  ' 


(25c) 


^  5  /I  V  I, 

R  a  _ 2_£. 

*^c®l  64 


__£ll 

2  2' 


K 

c  c 


(25d) 
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(25.) 


*cC' 


f 


S  ^  “pic 
64 


(25f) 


Substitution  of  Eqs.  (25a)  through  (25f)  into  Eqs.  (24a)  and  (24b) 
leads  to  the  following  expressions  for  to  first  order  in 
local  coordinate  systea  of  the  ring: 


ext 

6 


2ll 


Uo  T  p*cos4* 
4ii  ■'■c  R 

c 


3 


when  p'  <  And 


•  ext 
A0 


2 


%  J  p'cos^' 

4n  \  R 

c 


(26a) 


(26b) 


when  p'  ^  r^.  Let  (A,  a)  be  the  ring  position  in  the  local  cylindrical 
coordinates  of  the  toroidal  conductor  (cf.  Eqs.  (5a),  (5b)).  Then,  in  the 
relation  above,  R^,  p'  and  p'cosf'  are  replaced  by 


R  m  +  &  cos  a, 
c  o 


1/2 


p'  «  ^p^  +  -  2  p  A  cos  (♦  -  a)]  , 


p'  cos  4'  ■  P  CAS  4  -  A  cos  a. 


(27a) 


(27b) 

(27c) 
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Equations  (26a)  and  (26b)  are  useful  in  the  diffusion  problem  of  a 
current  ring  inside  a  toroidal  conducting  shell  vith  first  order  toroidal 
corrections.  This  vill  be  reported  elsewhere.  Here,  only  the  zero  order 
solution  to  the  diffusion  problem  is  considered.  In  this  case,  the  vector 
potential  of  the  current  ring,  in  the  absence  of  the  conductor,  is  equal  to 


.ext 

^eo 


(28a) 


inside  the  ring,  and 


.ext 

*9o 


2r  -^c 


(28b) 


outside  the  ring. 


.ext 

^00  -  “ 


o 

2r 


4r  ^c*'" 


Equation  (28b)  can  also  be  written  as: 

I  *  (;)  -  2  ; <♦-«)• 


(28c) 


This  expression  of  the  ring  vector  potential  is  used  in  the  application  of 
the  boundary  conditions  at  the  inner  surface  of  the  toroidal  conductor. 
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IV.  Honogeneous  Solution  of  the  Vector  Potential 

Inside  the  torus,  the  nost  general  solution  of  Bq.  (1)  is 

.in  .ext  .in 

where  is  the  homogeneous  solution  of  Eq.  (1).  The  exact  homogeneous 

solution  can  be  expressed  most  appropriately  in  toroidal  coordinates  with 

respect  to  the  toroidal  conductor  rather  than  the  current  ring.  For  that 

2  2  1/2 

purpose,  we  replace  Vt  M',  by  1^,  hj#  h^  -  Ir^  -«  1  *  in  Eqs. 

(10a),  (10b),  (11a),  (11b),  and  (11c).  When  hj  -  where  - 
£.n[(r^  +  a  +  h^)/(r^+  a  -  bj,)],  the  coordinates  (r,  2)  describe  a  circle 
which  coincides  with  the  inner  surface  of  the  toroidal  conductor.  In  terms 
of  the  toroidal  coordinates  (lij,  ^),  the  exact  honogeneous  solution  of  Eq. 
(1)  inside  the  torus  is 

m 

aJJJ  •  b^(coshr\|  -  cos^j)  ej,Q„  i  1/2  (coshn^) 

m>o 

♦  cosm^  +  sinm^j.  (30) 

For  small  aspect  ratio  a/r^,  and,  to  lowest  order  in  this  ratio,  we  have 
the  appproximate  relations: 

-r\.  -(n.  +  il.) 

bj  e  r^,  e  e  ®^2r^  and  e  «  pe  V2r^, 

where  (p,  ♦)  are  the  local  cylindrical  coordinates  with  respect  to  the 
toroidal  conductor  (cf.  Eqs.  (5a),  (5b)).  Moreover,  if  it  is  assumed  that 
each  of  the  coefficients  is  of  order  (r^/a)“,  then  it 

follows  from  Eqs.  (17b)  and  (30)  that  the  homogeneous  solution  of  Eq.  (1), 
to  xero  order  in  toroidal  corrections,  is  equal  to 
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(31) 


I 


■bO 


,1/2.  .  ^  3. 

n  r  (■  +  j) 

“T75 - 

(■  +  1) 


* 


Js) 

cosiD^  a_  siiiB# 

fli 


By  redefining  the  coefficients  1**  teres  of  the  zero  order 

coefficients  a  ,  a  ve  conclude  from  Eqs.  (29)  and  (31)  that  the 

most  general  solution  inside  the  torus  (inside  and  outside  the  current 
ring)  and  to  zero  order  in  toroidal  corrections  is 


DbI 

The  undefined  coefficients  a^,  a^^\  al®^  vill  be  determined  from  the 

OB  B 

boundary  conditions. 

The  homogeneous  zero  order  solution  satisfies  also  the  zero  order 
homogeneous  equation 


1  3_  ^^eho  1_  ^^^eho 

;  3p  ap  ^  1  3^2 


0, 


(33) 


which  follows  from  Eq.  (1)  by  expressing  it  in  the  local  coordinates  of  the 
toroidal  conductor  (cf.  Eq.  (6))  cuid  neglecting  the  terms  with  toroidal 
corrections. 

For  the  vector  potential  outside  the  torus,  we  define  the  toroidal 
coordinates  (h^,  ^^)  in  a  similar  fashion,  i.e.,  we  replace  ,  h'  and  b^ 
by  b^  -  [r^^  -  in  Eqs.  (10a),  (10b),  (11a),  (11b)  and 

(11c).  When  -  tlQj.,  where  +  b  +  bjj)/(rQ  +  b  -  b^)l,  the 

coordinates  (r,z)  describe  a  circle  which  coincides  with  the  outer  surface 
of  the  toroidal  conductor.  In  terms  of  the  toroidal  coordinates  (h^f  ^g)* 
the  exact  solution  of  Eq.  (3)  outside  the  torus  is 
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‘‘T  ’  K  Y.  -  1/  2 

m*o 

As  before,  for  small  aspect  ratio  b/r^,  ve  have  the  approximate  relations: 

-ri  'Irt  4  A 

b^  «  r^,  e  «  b/2r  and  e  °  °  «  2r  e^Vp.  If  It  Is  assumed  that 

0  0  o  o 

each  of  the  coefficients  b  b^S_^®^  is  of  order  (b/r  )",  then  it 

O  B  O  H  O 

follows  from  Eqs.  (17a)  and  (34)  that  the  solution  of  Bq.  (3),  to  zero 
order  in  toroidal  corrections,  is  equal  to 

m-l 

Here,  the  coefficients  have  been  redefined  in  terms  of  the  zero 

B  B 

order  coefficients  b^,  bl®\  bf®\  These  undefined  coefficients  will  be 

OB  B 

determined  from  the  boundary  conditions. 

It  is  apparent  that  Eq.  (35)  satisfies  the  zero  order  homogeneous 
equation  (33),  but  the  solution  of  Eq.  (33)  does  not  provide  all  the 
information  included  in  Eq.  (35).  The  most  general  solution  of  Eq.  (33) 
which  is  independent  of  the  toroidal  angle  4  is  equal  to  •«-  (.np  where 
and  are  arbitrary  constants.  But  Eq.  (35)  indicates  that  these  two 
constants  are  related  to  each  other  and  their  dependence  on  each  other  is 
established  only  by  solving  the  problem  in  toroidal  geometry  rather  than 
making  some  ad  hoc  assumption.  For  example,  if  we  assumed  that,  at 
infinite  time,  the  vector  potential  outside  the  conductor  is  equal  to  that 
of  the  current  ring  in  the  absence  of  the  conductor,  ve  would  probably 
obtain  the  correct  relationship  between  and  C^,  but  this  assumption 
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vould  be  iaposed  on  the  solution  of  the  diffusion  problem  rather  than 
coning  out  naturally  as  a  result  from  the  solution. 
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V.  -Initial  Conditions  of  the  Vector  Potential 

When  t  and  a  are  tine  dependent  quantities,  the  coefficients  a^, 

Eqs.  (32)  and  (35)  are  also  tine 

H  H  O  D  n 

dependent.  Since  outside  the  torus  the  vector  potential  is  zero  at  t  >  0, 
h^(t),  b£^^(t)  and  bl^^(t)  are  also  zero  at  t  ■  0.  In  addition,  the  vector 
potential  is  zero  at  t  -  0  inside  the  conductor.  Fron  the  continuity  of 
the  vector  potential  at  the  inner  surface  of  the  conductor  and  fron  Eq. 
(32),  it  follows  that  the  coefficients  a  (t),  a^^\t),  ai^^(t)  are  not  zero 
at  t  «  0.  Since  the  image  fields  constitute  a  zero  order  homogeneous 
solution  inside  the  torus,  it  is  convenient  to  redefine  the  as  yet 
undetermined  coefficients 

f  C  ^  /  e  \ 

ajl^  (t),  ajl^  ^(t)  by  subtracting  the  image  solution  from  them,  so  that  they 
are  zero  at  t  •  0.  As  to  the  coefficient  a^(t),  in  order  that  it  becomes 
zero  at  t  -  0,  it  is  convenient  to  redefine  it  by  replacing  it  with  a^(t)  - 
(y^/2ii)I^(t)  (tn  8r^/a  -  2).  Then,  Eq.  (32)  should  be  replaced  by 

_  ^  A^(t)  -  2p  A  (t)  cos  (♦  -  a(t)) 

■►fen  fl  -  2  cos  (♦  -  «(t))] 

a  a 

m 

+  a^(t)  +  ^  (l)  cosm^  a^®\t)  sinm^], 

m«l 
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inside  the  ring,  and 


aJJ  ^  I^(t)  [j  tn  I  (36b) 

-2S^c.,(6  -  a(t))) 

+  -  2  cos  (♦  -  a(t))]] 

ft  ft 

m 

+  a^(t)  +  ^  (j)  cosa^  +  a^®^(t)  sina#], 

■■I 

outside  the  current  ring  but  inside  the  torus.  The  vector  potential  froa 
the  iaage  has  a  logarlthaic  singularity  at  the  iaage  position  (a  /A(t), 
«(t)),  which  lies  outside  the  inner  surface  of  the  conductor.  Therefore, 
inside  the  torus,  it  is  a  zero  order  hoaogeneous  solution  of  Eq.  (33). 
When  d(t)/p  <  1,  £q.  (36b)  is  equivalent  to 

■  n  V*)  [‘"  I 

-  i  h  m"  m"  -  fen  -  ■  <♦  -  •<«»] 

m-1 

+  a^(t)  +  ^(j^)  cosn4  +  a^®^(t)  sina+j,  (37) 

m«l 

where  a  (t),  a^‘^^(t),  a5®^(t)  are  zero  at  t  «  0.  The  zero  initialization 

O  IB  D 

of  these  coefficients  will  lead  to  siaple  expressions  when  the  boundary 
conditions  will  be  applied  in  the  next  Section. 
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VI.  -Diffusion  Fields  Inside  and  Outside  the  Toroidal  Conductor 

In  order  to  compute  the  fields*  the  zero  order  diffusion  equation  aust 
be  solved  inside  the  conductor  and  the  boundary  conditions  be  applied  on 
its  inner  and  outer  surface.  The  zero  order  diffusion  equation  is  obtained 
froa  Eq.  (6)  by  omitting  the  toroidal  corrections*  namely, 


1  9 
P  ap  ^ 


dA 


con 

6o 


dp 


con 


1  ^^eo 

7 


dA 


»'o" 


con 

60 


at 


0 


In  general,  let 


(38) 


f<p) 


J  f(t)  e'P*  dt. 


(39) 


be  the  Laplace  transform  of  f(t).  Then,  in  the  Laplace  transform  domain, 
Bq.  (38)  becomes 


con 


1  a 

p  ap  ^ 


»^6o 

fT" 


1 

7 


39^ 


-Up 

**o" 


p  Ag"  -  0, 


(40) 


where  A^"  (p,  9,  t)  was  assumed  to  be  zero  at  t  ■  0.  The  most  general 
solution  of  Eq.  (40)  is 


A^"  (P.t.p)  -  c^(p)  VXp) 


m-1 


♦  d^(p)  Ko(Xp)  + 


^  K„(Xp)  cosat  +  d^®^(p)  sinat],  (41) 

m«l 
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where 


X  .  [l.„.  p]  .  (42) 

and  I_(x),  K_(x)  are  the  modified  Bessel  functions  of  order  a. 

ID  01 

The  boundary  conditions,  to  zero  order  in  toroidal  corrections,  follow 
from  Eqs.  (7a),  (7b),  (8c),  (8d),  (9a),  (9b)  and  (39).  They  are  given  by 


Ai"  (a,*,p)  -  A^"(a,^,p), 


(43a) 


aA^(p>».p) 


aA|g^(p,4,p) 


(43b) 


and  the  same  two  relations  at  p«b.  Defining  and  F^®^(p)  to  be  the 

Laplace  transforms  of  (A(t)/a)'"  cosm  a(t)  and  (d(t)/a)"  sinm  a(t), 
respectively,  the  boundary  conditions  at  p-a  and  p>b,  using  Eqs.  (35),  (37) 
and  (41),  lead  to  the  following  algebraic  system  of  equations: 


ao<P)  -  c^(p)Io(Xa)  ♦  d^(p)K^(Xa), 


(44a) 


[in  ^  -  2]b^(p)  -  c^(p)I^(Xb)  +  dj^(p)Kjj(Xb),  (44c) 

-  b^(p)  -  Ab[c^(p)i;  (Xb)  +  d^(P)K;  (Xb)],  (44d) 

when  m-0,  and 

‘i*’  (P)  -  *  4i‘^P)K„(Xa),  (45.) 
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(45b) 

-  n  ^a\?)  -  ^fi^^P)i;  (W»  +  d^^\p)K;  (W))],  (45d) 

when  m  >  1,  2,  ...,  and  i  ■  c,  s.  Also,  Il(x),  Rl(x)  are  the  derivatives 

B  B 

of  Ig,(x),  Kg|(x).  The  solution  of  the  first  algebraic  system  of  Eqs.  (44a) 
-  (44d)  is 


* 

Pq  * 

La) 

%<P' 

-  H  ^c<P> 

(a> 

b^P) 

•  H  ic(P> 

1 

foCX) 

Co(P) 

•  Vp)  [( 

tn  -gS  -  2jXbKj(Xb)  -  K^(Xb)J, 

8r  »  ^ 

do<P) 

-  b^(P)  [( 

tn  -g2  -  2jXb  Ij(Xb)  +  Io(^)J» 

where 


(46a) 

(46b) 

(46c) 

(46d) 


fp(X)  -  (tn  -  2)XaXb  [Kj(Xa)Ij(Xb)  -  Ij(Xa)Kj(Xb)] 

+  Xa[K^(Xa)I^(Xb)  +  Ij(Xa)Kjj(Xb)], 

8ir 

go(X)  -  (tn  -gS  -  2]xb[K^(Xa)Ij(Xb)  +  I^(Xa)Kj(Xb)] 


(47a) 


(47b) 


+  K^(>a)I^(Xb)  -  I^(Xa)K^(Xb). 

Siailarly,  the  solution  of  the  second  algebraic  system  of  Eqs.  (45a)  - 
(45d)  is 


(P) 


r.(X)’ 


(48a) 


f,<X) 


ci‘’(P)  -  b<*><p)  Xb  K^i<Xb), 
-  bi‘>(P)  Xb  I^,<Xb). 


where 


(48b) 

(48c) 

(48d) 


r,(X)  .  -  ^  [l,*,<»«)*:^l(Xb)  -  <!„i(X*)I^l(Xb)],  (49.) 

g^(X)  -  Xb  [l,(Xa)K^_j(Xb)  +  K^(Xa)I^_j(Xb)],  (49b) 

and  m  -  1,  2,  3,  i  -  c,  s. 

The  inverse  Laplace  transforms  of  the  coefficients  given  above  are 

determined  from  the  sum  of  the  residues  at  the  poles  of  these  coefficients. 

The  poles  are  computed  at  the  zeroes  of  ^g|(\|]^)  *  0,  where  a  -  0,  1,  2, 

...»  k>0p  1,  2,  ...  Since  all  the  zeroes  occur  for  imaginary  values  of 

X,  we  define  the  real  quantities  by  means  of  the  relation  ^mk  -  i»mk- 

2 

Then  a  pole  occurs  at  p^  *  inverse  Laplace 

transform  h  (t)  of  g„(X)/?  (X)  is 
m  m  m 
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Also, 


(52c) 


(53a) 


(53b) 


V  “  -  «^f;(«^) 


where  a>0,  1,  2,  ....  k> 


g.(a)  and  the  derivative  f' 

D 


(53c) 

>  0,  1,  2,  ...  Finally,  the  functions  £  (a) 

B 

_(a)  of  f_(«)  are  as  follovst 

B  B 


f 


*o<“>  •  <“*> 

8r 

*.<«)  ■  r.<i«)  ■  I 

(54b) 

g„<.)  ■  g„(i.>  -  f  -  ’'o<’‘o>-'o<’'l>]  <“•> 

8]r 

«„(«)  -  «„(!«)  -  f  *l[V’'o>Vl<*l>  -  V’'o>Vl<’'l>]’  <55b) 
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(56a) 


of' (a)  ■  iaf'(ia) 

■  ‘o<«>  *  ?  -  ’'o<-o>-’o<-l>] 

Sir 

Sit 

-  H(‘"  -r  -  2)>'?*l]’‘o['>X<-o>’'o<-‘l>  -  ’x<*.>-’o«-X>] 

8r 

~  2  (*•"  "b  ^]*o*l 
ofjjjCa)  ■  ia  f;(i«)  (56b) 

X^x 

-  f  -V  -  V’'o>Vi<*i)] 

X  x^ 

-  f  -!^  [Vx<='o>’'.<’'x>  -  Vx<»o>V*x>]- 

where  x  -  oa,  x.  ■  ob,  a  •  1,  2,  ...  and  J.(x),  T_(x)  are  the  Bessel 
functions  of  order  m.  Notice  that  are  the  zeroes  of  ^j,(“ijjf^)  “  0» 

®*0j  If  2f  «..f  lc*0|  If  2f  ... 

From  Eqs.  (36a)f  (36b)f  (35)f  (Sla)-(51d)f  we  conclude  that  the  zero 
order  vector  potential  inside  the  ring  is  equal  to 


c 
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+  tn  1  +  -  2  cos  (♦  -  a(t))jj 

•  «B  OB 

♦  Z*ok"ok<*>  *  Z  Z  V©'  [“‘iitt)'”"*  * 

k-O  1^  k«0 

Outside  the  ring  and  inside  the  torus  <i.e,  p  ^  a),  it  is  equal  to 

^  (57b) 

-  tn  (l  +(^)^  -  2  ^  cos  (♦  -  a(t))) 

+  tn(l  +  -  2  cos  (i  -  e(t)))] 

A  ft 

•  mm 

*  l*<,k»0k<'>  ♦  I  Zv©* 

mZi  kmO 

and  outside  the  torus  (i.e«  p  ^  b)  it  is  equal  to 

8r  * 

-  2)yBok»„k<')  <“> 


*  Z  Z‘-k©”  “»♦*  <*>  »‘H- 

a-l  k-0 

From  Bqs.  (52a),  (52b)  and  (52c)  it  is  easy  to  show  that  the  time  dependent 
coefficients  U^i^(t),  U^^(t)f  satisfy  the  first  order  differential 

equations 


«ok<‘>  ‘  ^  “ok<')  -  ^  Ji 


(59a) 
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(59b) 


U<j)(t)  *  ^  oM(t)  .  ^  ^  Vt)  (^)‘  com-rtt) 

VI<^>(t)  *  U<J>(t)  .  ^  I^<t)  (^]"  sin..<t)  (590 

with  initial  conditions  U^j^(o)  ■  ■  ■  1» 

2,  ...,k  *  0,  1,  2,  ...  These  differential  equations  are  very  useful  when 
the  current  ring  moves  and  its  equations  of  notion  depend  on  the  diffusion 
fields,  l.e.,  vhen  the  ring  dynamics  is  coupled  to  the  diffusion  fields. 
Then,  the  state  vector  of  the  system  consists  not  only  of  the  position  and 
velocity  of  the  ring,  but  also  of  the  diffusion  coefficients  Uo^(t), 
U^^(t),  U^\t),  for  n  -  1,  2,  ...»  k  -  0,  1,  2,  ...,  and  its  tine 
derivative  is  determined  by  the  ring  dynamics,  as  veil  as  by  Eqs.  (59a)- 
(59c).  Notice  that  it  is  much  easier  to  solve  in  the  computer  a  set  of 
coupled  first  order  differential  equations  rather  than  a  set  of  coupled 
first  order  differential  equations  and  the  convolution  Integrals  given  by 
Eqs.  t52a)-(52c). 

Next,  four  exact  identities  will  be  established  for  the  time 
independent  coefficients  For  a  motionless,  step  function  ring 

current,  i.e.,  when  I^.(t)  -  0(t),  Eqs.  (51a)-(51d)  and  Eqs.  (52a)-(52c) 

give 


a„(t) 


bo(t) 


(60a) 


(60b) 
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*  i  ^  I.g)*  ^  .  .■"^•>'),<«,c) 


*  1  bl"(t)  .  :^I.|-|  ■ 


i»(t)  -  j|  .*•“ 


-t/T 


(l  -  »  (60d) 


II  -  . 


The  Laplace  transform  of  ^^(t)  is  I^(p)  >  ve  have  from  Eqs. 

(46a),  (46b),  (48a),  (48b): 


P.  ,  g.(X) 


«  -2  T  1  ^ 

®0'P)  “  9n  » 


2n  o  p  f 


(61a) 


b  (p)  .  -O  I  i  _1 - 


(61b) 


•m  <P>  ♦  <P>-H^ctJ  *  Pil^’ 


(61c) 


*  t  i.<P)  ■  jS  1,  (?)  I 


(61d) 


From  the  vell-knovn  theorem  of  Laplace  transforms  it  follows  that  ^^§f(t) 
pf(p).  Eqs.  (47a),  (47b),  (49a),  (49b),  give 


kiS  ^o^^  )  - 


(62a) 


8r. 


Jlig  go(X)  -  tn  —  -  2, 


(62b) 


iig  f„(X)  -mg). 


(62c) 
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Vlth' the  help  of  . these  identities  and  Eqs.  (35),  (36a)  and  (36b),  it  is 
easy  to  show  that,  as  t  -»  •,  the  zero  order  vector  potential  inside  and 
outside  the  toroidal  conductor  (as  veil  as  Inside  the  ring)  becoees  equal 
to  that  of  the  current  ring  in  the  absence  of  the  conductor,  i.e.,  it 
becones  equal  to  ,  due  to  the  diffusion  process.  This  conclusion 
demonstrates  the  importance  of  the  identities  (63a)-(63d). 

The  magnetic  and  electric  fields  inside  and  outside  the  torus  can  be 
computed  from  Eqs.  (57a),  (57b)  and  (58).  Thus,  the  self-magnetic  and 
self-electric  fields,  i.e.,  the  fields  of  the  ring  at  its  centroid  p  > 
h(t),  ^  -  a(t),  are  equal  to 


■  1 1  !  V 


(65a) 


m«l  k-o 


*  [-  U^\t)  cosmoCt)  .  U^>(t)  sliw.(t)  ], 


-  V«)  ^ 


(65b) 


1 1  iv  (^r 


m«l  k>iO 


*  [^mk^^^^  cosma(t)  +  U^J^(t)  sinma(t)  ], 

■  -  'll  [-  f:  *  I  -n  (l  -  i^f)] 


*  H  V') 


Xc(t)Xc(t)  2^(t)Z^(t) 


1  -  (^)- 
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(65c) 


-  cosma(t)  -  U^^(t)  sln«a(t)jl, 

where  x^(t)  -  h(t)cosa(t)  and  z^(t)  -  A(t)  sina(t). 

self 

The  last  three  relations,  and  in  particular  the  component  were 

useful  in  providing  an  analytical  model  for  the  beam  trapping  that  occurs 
after  injection  in  the  modified  betatron  accelerator. 

A  typical  example  of  the  effect  that  the  diffusion  process  has  on  the 
ring  dynamics,  immediately  after  injection  is  shown  in  Fig.  2.  The 
numerical  integration  of  the  ring  equations  of  notion  coupled  with  the 
diffusion  fields  was  done  for  the  parameters  listed  in  Table  I.  Figure 
2(a)  shows  the  projection  of  the  centroid  orbit  on  the  r-z  plane  that  moves 
with  the  same  toroidal  angular  velocity  as  the  ring  centroid.  There  is  a 
slow  (bounce)  motion  and,  due  to  the  presence  of  the  stellarator  windings 
(i.e.,  strong  focusing),  the..e  is  also  an  intermediate  motion.  Both  of 
these  modes  are  indicated  in  Fig.  2(a).  Since  there  are  six  field  periods 
of  the  stellarator  field  in  the  range  o  <  6  <  2n,  the  electrons  perform  six 
oscillations  during  one  revolution  around  the  major  axis.  To  take  into 
account  the  intermediate  motion  that  has  been  neglected  in  the  diffusion 
model  presented  in  this  paper,  the  resistivity  in  the  code  is  computed 
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using  the  skin  depth  that  corresponds  to  the  interaediate  frequency  and  not 
the  actual  thickness  of  the  %rall.  The  dots  in  Fig.  2(b)  show  the  positions 

o 

the  beu  crosses  the  r-z  plane  at  6  •  0  .  The  tine  difference  between  two 
dots  is  equal  to  the  period  around  the  aajor  axis,  i.e.,  -  27  nsec,  and 
therefore  the  speed  of  the  ring  on  the  r-z  plane  can  be  inferred  from  the 
relative  position  of  the  dots.  Fig.  2(c)  provides  the  relativistic  factor 
Y  vs.  tine  and  the  reduction  in  r  is  obvious  due  to  the  energy  lost  on  the 
resistive  wall  and  to  establish  the  electromagnetic  field  outside  the 
torus.  Another  example  is  given  in  reference  13,  which  refers  to  the  beam 
trapping  in  the  modified  betatron  accelerator  and  is  in  good  agreement  with 
the  experimental  results. 
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VII..  Approxiaate  Results  for  a  Thin  Conducting  Vail 

The  results  presented  in  the  previous  Sections  are  approxiaate  in  the 
sense  that  they  include  only  zero  order  toroidal  corrections,  i.e.,  are 
valid  only  for  soiall  aspect  ratio  vessels.  Othervise,  they  are  exact.  In 
this  Section,  the  additional  assuaption  is  aade  that  the  conducting  vail  is 
thin.  This  assuaption  allows  us  to  coapute  approxiaate  expressions  of  the 
zeroes  of  ■  0  (cf.  Eqs.  (54a),  (54b))  and  of  the  vector 

potential  and  the  fields. 

When  the  conducting  vail  is  thin,  i.e.,  (b-a)/a  «  1,  both  x^  ■  ea  and 
Xj^  ■  ed),  where  a  is  a  zero  of  f|||(a)  -  0,  are  very  large  nuabers  and  the 
asyaptotic  expansions  of  the  Bessel  functions  can  be  used  in  Eqs.  (54a), 
(54b).  This  is  valid  only  up  to  soae  aaxiaua  value  of  a.  Keeping  terns  up 
to  order  1/z,  the  asymptotic  expansions  of  ^^(z),  T^(z)  are: 

■'»<')  -  JF  *»]•  <“•> 


-JFt 


.1  . 

sin  "  co®  K 

n  oz  r 


(66b) 


where  «  z  -  (n  1/2)  ii/2.  Ve  substitute  these  expansions  into  Eqs. 
(54a),  (54b).  Then,  the  zeroes  of  determined  from 

tan(Xj  -  x^)  .  - — jjl - — .  (67.) 

(tn  -jS  -  2)x„ 


while  the  zeroes  of  £g,(«p||^)  ■  0  are  determined  from 


tan(Xj  -  x^) 


(67b) 
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A  more  accurate  expression  for  the  zeroes  of  ■  0,  correct  to  order 

2 

(b/x^)  ,  is  derived  in  the  Appendix,  and  is  given  by 


tan(Xj  -  x^) 


1  - 


r 


(67c) 


If  x^  »  1,  ve  see  from  Eq.  (67a)  that  |xj^  ~ 
tan(Xj^  -  x^)  ■  x^  -  x^  a  (b  -  a)a,  and  one  of  the  zeroes  is 


oo 


[(tn  -  2)a(b  - 


a) 


1/2 


(68a) 


while  the  others  are  given  by 


*ok 


kn 

b  -  a’ 


(68b) 


where  k  -  1,  2,  ...  The  snail  additive  correction  tern  l/[kii(tn(8r^/b) 

-  2)a]  has  been  omitted  in  Eq.  (68b).  The  zeroes  of  Eq.  (67b)  can  be 
obtained  in  a  similar  fashion,  except  when  m  is  as  large  or  larger  than  x^. 
Let  m  -  1,  2,  ...,  M,  where  N  -  Int(a/4(b  -  a)],  and  Int(x)  is  the  integral 
part  of  X.  Then  one  set  of  zeroes  of  Eq.  (67b)  is  approximately  given  by 


V  •  tr^rrr]  • 

_2 

with  an  error  of  only  a  few  percent  when  (b  -  a) /a  <  10  .  For  m  as 
specified  above  and  k  >  1,  2,  ...,  K,  where  K  >  Int[a/n(b  -  a)],  the  rest 
of  the  zeroes  are  given  by 


*mk 


kn 

b  -  a* 


(69b) 
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The  small  correction  term  [2m  +  (kn)  /2]/(kna)  has  been  omitted  in  Eq. 

(69b).  For  values  of  a  and  k  larger  than  M  and  K,  respectively,  the  zeroes 

should  be  computed  numerically  from  Eq.  (54b).  Notice  that  the  presence  of 

the  terms  (p/a)"*  and  (b/p)"*  in  the  series  expansions  of  the  vector 

potential  and  the  fields  Indicates  that  the  large  values  of  m  become 

important  when  these  quantities  are  computed  close  to  the  conducting  vail, 

where  p/a  and  b/p  become  almost  equal  to  1  and  more  m-terms  must  be 

included  in  the  sums  to  converge  vithin  a  prescribed  accuracy.  An  estimate 

of  the  minor  radius  p^^  vithin  which  the  vector  potential  and  the  fields  are 

sufficiently  accurate  is  determined  by  p^/a  >  [(b  -  a)/a]  .  Since 

(b  -  a)/a  «  1,  all  the  terms  associated  with  (p/a)"*,  for  m  ■  M  1,  N  2, 

...  in  the  series  expansions  of  the  vector  potential  and  the  fields,  have  a 

negligible  contribution,  provided  p  <  p^.  A  similar  argument  can  be  made 

for  the  vector  potential  and  the  fields  outside  the  torus.  Their  accuracy 

1/M 

is  vithin  a  fev  percent  when  p  ^  p^,  where  p^/a  -  [a/(b  -  a)]  .  Vithin 

the  distances  d^  >  a  -  p^  and  ■  Pq  ~  ^  inner  and  outer 

conducting  vails  the  zeroes  (and,  therefore,  the  vector  potential  and  the 
fields)  cannot  be  computed  analytically  in  terms  of  a  simple  expression. 

In  this  case  they  should  be  computed  numerically  from  Eq.  (54b)  and  then 
use  the  analytic  expressions  for  the  vector  potential  and  the  fields.  As 
an  example,  when  (b  -  a)/a  «  10~^,  then  M  >  250,  K  «  318,  d^/a  «  0.027, 
d^/a  *  0.028,  but  when  (b  -  a)/a  »  10”^,  then  M  »  25,  K  «  31,  d^/a  =  0.17, 
d^/a  ■■  0.20.  In  the  following,  the  various  quantities  will  be  computed  to 
order  (b  -  a) /a,  with  the  understanding  that  they  are  not  accurate  close  to 
the  conducting  vail.  But  in  the  limit  when  the  ratio  (b  -  a)/a  tends  to 
zero  (but  <r(b  -  a)/a  remains  finite),  i.e.,  when  the  toroidal  conductor 
becomes  a  toroidal  conducting  shell,  the  distances  d.,  d^  are  zero  and  the 
results  become  exact,  to  zero  order  in  toroidal  corrections,  everywhere 
inside  as  well  as  outside  and  in  the  vicinity  of  the  torus. 


-..Under  the  assuaptions  stated  above,  the  tine  constants  in  Eq. 
(53a)  becoae 


oo 


V- 


2  ^)* 

(70a) 

(70b) 

**1) 

a  * 

(70c) 

*'PD 

k2  ’ 

(70d) 

vhere 


U  o(b  -  a)a 

“  2  ’  (71a) 

M  e(b  -  a)^ 

j - ,  (71b) 

n 

and  a  «  1,  2,  k  >  1,  2,  ...  Therefore,  there  are  three  characteristic 
tiae  constants  associated  vith  the  diffusion  process.  The  "loop"  diffusion 
tiae  is  the  slowest  and  deteraines  the  speed  with  which  the  external 
field  of  the  ring  diffuses  into  the  hole  of  the  doughnut.  This  tiae  is 
present  because  of  the  toroidal  geoaetry  of  the  conductor.  The  "cylinder" 
diffusion  time  and  the  "fast"  diffusion  tine  are  associated  with  the 

diffusion  process  in  a  cylinder  and  deternine  the  speed  with  which  the 
field  of  the  ring  penetrates  the  conducting  wall.  Notice  that,  in  the 
limit  of  a  toroidal  conducting  shell,  any  terms  associated  with  the  fast 
diffusion  tine  Tpi^  diffuse  instantaneously  at  t  «  0  outside  the  torus. 
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This,  explains  the  origin  of  the  electric  field  that  is  iaaediately 
established  at  t  >  0  outside  the  torus  for  the  shell  eodel  (cf.  stateemt 
after  Eq.  (90)  in  the  next  section).  On  the  other  hand,  we  know  free  Eq. 
(58),  that  the  magnetic  field  does  not  diffuse  instantaneously  at  t  «  0 
outside  the  torus. 

Under  the  same  assumptions  stated  above  and  to  lowest  order  in 
(b  -  a)/a,  we  have  the  following  approximate  relations 


A  e  tn  —  -  2  -  y  — - — , 
oo  a  j  a 


A  «  2  b  -  a 

(kn)^  ® 


A  «  I  2 

mo  m  3  a 


A  e  4  b  -  a 

*nik  /•  \2  a  ’ 
(kn) 


oo 


b  -  a 


B  u  * 
ok 


2  (-1)' 


b  -  a 


8r.  .  2  ® 


(en  -gS  -  2)(kn)^ 


(72a) 

(72b) 

(72c) 

(72d) 

(73a) 

(73b) 


B 

mo  m  3  a 


(73c) 


B  ^ 

(kii)2  » 


(73d) 


for  m  -  1,  2,  k  «  1,  2,  ...  The  relations  (73b),  (73d)  were  derived 

directly  from  Eqs.  (56a),  (56b),  while  the  relations  (72b),  (72d)  were 
derived  with  the  help  of  Eqs.  (73b),  (73d)  and  the  relations  * 


(-1)*^  (t-n  8r^/a  -  2),  ®  (-1)*^»  m  -  1,  2,  . . . ,  k  «  1,  2,  . . . 
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Finally,  were  derived  using  the  identities  (63a)  - 

(63d). 

As  mentioned  above,  the  terns  associated  with  the  fast  diffusion  tine 
Tpjj  vary  in  tine  on  a  nuch  faster  tine  scale  than  the  terns  associated  with 
the  tines  and  If  the  ring  current  l^(t)  and  its  position  (A(t), 
a(t))  vary  slowly  within  a  few  e-folds  of  Xpj^,  then  the  part  of  the  vector 
potential  (or  the  fields)  which  is  associated  with  Xp^  can  be  simplified 
considerably. 

First,  let  us  consider  the  self-nagnetic  field.  Substitution  of  Eqs. 
(70b),  (70d)  and  Eqs.  (72a)  -  (72d)  into  Bqs.  (65a),  (65b)  leads  to  the 
relations 
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2m  b  -  aj 


(74a) 


*  [”^no^^^^  cosna(t)  +  sinma(t)j 
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*  cosma(t)  +  U^®^(t)  sinmo(t)] 
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where 


t'  2 

f *-''>  •'  ^  ««'• 


h”“(t,  t')  .  Hl^(t')  44^ 

[1  ,  eo.  (.(,)  -  .(.'))  -  2 

[1  .  (4itl|Ull.)"  .  2  Miipi  „.  (^.>  .  ^,,,)]" 


(75a) 


hf®^^t,  t')  -  ^  1^(1') 


(75b) 


[1  ,  [^Lnmilf]  sin  («(t)  .  a(t')) 

ft 

[1  .  (MIl|Ua)"  .  j  iiSJfiD.  CO,  (.^.)  .  ^,.))]' 


If  I^(t  -  t')»  A(t  -  t'),  a(t  -  t')  vary  only  slightly  as  t'  varies  within 
a  few  e-folds  of  Tp^,  they  can  be  replaced  by  I^(t),  h(t)  and  a(t),  except 
in  the  sine  in  Eq.  (75b)  we  should  set  a(t)  -  «(t  -  t')  »  a'(t)t'  to  get 
the  lowest  order  contribution.  Here  a'(t)  is  the  derivative  of  a(t). 
Equations  (74a)  and  (74b)  then  becoae 


0).ix  (*-T*^)(^r 

iii>i 

*  f-  cosna(t)  +  ui^^(t)  sinna(t)] 
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(76b) 


2 

After  a  few  e-folds  of  Xp^,  the  sub  over  k  becoaes  equal  to  n  /6  and  there 
is  a  residual  contribution  in  the  self-aagnetic  field  froa  the  fast 
diffusing  terns.  Since  (b  -  a)/a  «  It  Eqs.  (76a) ,  (76b)  indicate  that 
this  contribution  is  snail  unless  the  ring  is  close  to  the  conducting  wall. 
But  in  that  region,  these  relations  are  no  longer  valid  and,  therefore, 
they  provide  only  a  hint  as  to  the  significance  of  the  fast  diffusing  terns 
when  the  ring  is  close  to  the  conducting  wall. 
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■Froa  Bqs.  (65c),  (70b)  (70d)  and  Eqs.  (72a)  -  (72d),  and  to  lovast 
order  in  (b  -  a)/a,  the  electric  field  at  the  ring  centroid  is  given  by 
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VhMi. A(t),  a(t)  vary  slovly,  in  a  siailar  fashion  as  for  tha  salf- 
■agnatic  fiald,  va  obtain  tha  siaplifiad  aquation  for  tha  salf-alactric 
fiald 


psalf 

®eo 


<t) 


-  2i 


1 

I* 


In 


♦ 


**0 

Jr 


Ie<t) 


(^)  Ko^<^>  cos««(t)  ♦  uij^t)sin«a(t)] 

■>1 


»*o 

-  Tj, 


©3(0, 


vhara  6^ 


(s,  q)  is  tha  thata  function  of  order  3,^^  i.a., 

©3(z«  q)  -  1  +  2  y  q  cos  2  k  z. 

1^ 


(79) 


(80) 


Notice  that  0^  (0,  q)  can  be  expressed  in  terns  of  the  conplete  elliptic 

1/7 

integral  of  the  first  kind  K(b),  i.e.,  62(0,  q)  >  [2  K(B)/n]  .  Here,  ■ 

is  related  with  the  none  q  by  the  relation  q  -  exp[-n  R(l-n)/R(B)]  and  when 
q  tends  to  1,  then  ■  tends  also  to  1.  When  t/Tpj^  «  1,  the  none  q  is  very 
close  to  1,  and  in  this  case  R  (!-■)  s  R(o)  «  n/2,  so  that 

[2R(m)/Rj^^^  «  [R/tn(l/q)]^^^.  Therefore,  83(0,  e  ^^^®)  «  [n/(t/Tpj,) 

1/2 

i.e.,  the  self-electric  field  is  proportional  to  I^(t)[TpQ/t]  vhen  t/Xpp 
«  1.  If  the  ring  current  is  a  step  function  of  tine,  i.e.,  l^(t)  > 

1^8(1),  then  the  self-electric  field  is  infinite  at  t  «  0  (actually,  it  is 
infinite  everywhere  inside  the  torus).  This  result  is  not  surprising  if  we 
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tak*.  into  account  that  the  vector  potential  rises  in  tine  instantaneously 
at  t  -  0,  because  the  ring  current  does  so.  When  t/Xpp  »  1,  ve  have  83(0, 

e  )  B  1  and  the  fast  diffusing  part  in  provides  a  residual 

contribution. 

Finally,  consider  the  electric  field  outside  the  torus,  i.e.,  vhen 
P  ^  b.  From  Eqs.  (4c),  (58),  (59a)-(59c),  (70d)  and  Eqs.  (72a)-(72d)  we 
see  that,  to  lowest  order  in  (b-a)/a,  it  is  given  by 
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where 


hQ(p,4,t) 


bd(t) 

pa 


(82) 
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*  (♦  -  *<^>]  - 

When  I^(t),  A(t)»  a(t)  vary  alovly,  Bq.  (81)  siaplifias  to 


.out 

®eo 


(pf4»t) 
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■  "oo(^> 
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(83) 
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^  Z  S)"  •*“♦) 
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vherc  6^(z,  q)  is  the  theta  function  of  order  4,^®  i.e., 


e 


^(z,  q)  -  1  +  2  ^  (-!)•"  q*^^  cos  2kz. 


(84) 


Notice  that  0^(0,  q)  -  I(l-a)^^2  2  K(a)/nJ^^2,  when  t/Xp^  «  1,  the  none  q 
is  very  close  to  1  and  in  this  case,  R(m)  e  (1/2)  en[16/(l-m) ]  and  2  R(b)/ii 
s  n/(.n(l/q).  Therefore,  ve  have 


It 


t/x, 


It  2(t/Xpjj) 

e 


FD 


1/2 


(85) 


i.e.,  vhen  the  ring  current  is  a  step  function  of  tine,  the  electric  field 
is  zero  at  t  -  0  outside  the  torus.  On  the  other  hand,  vhen 
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»  1,  then  0,e  «  1  and  the  electric  field  reduces  to 

that  computed  from  the  shell  model.  It  appears,  therefore,  that  the  fast 

diffusing  part  in  contributes  the  exact  amount  needed  for  eJ;'  to 

vanish  at  t  >  0,  but  very  quickly  increases  to  the  value  predicted  by 

the  shell  model  vhen  the  conducting  vail  is  thin. 

To  calculate  the  wall  current  density  J^^^(p,4,t)  -  t),  the 

electric  field  inside  the  conductor  is  needed.  From  the  continuity  of  the 

electric  field  at  the  inner  and  outer  surfaces  of  the  conductor  and  vhen 

t/Tpij  «  1,  E^"  varies  from  a  very  large  value  at  the  inner  surface  to  a 

very  small  value  at  the  outer  surface.  However,  vhen  t/Tp^^  »  1,  but  t/x^ 

«  1,  and  in  the  special  case  of  a  thin  conducting  vail,  it  is  easy  to  show 

from  Eqs.  (S7b),  <72a),  (72c)  and  (83)  that  the  electric  fields  at  the 

inner  and  outer  surface  are  approximately  equal  to  each  other.  In  the 

extreme  case  of  the  shell  model,  they  become  exactly  equal  to  each  other, 

wftXX 

and  the  surface  vail  current  density  is  equal  to  (4,t)  >  o(b 
a)E^*(a,4,t),  where  E^^(a,+,t)  is  given  by  Eq.  (83).  The  surface  vail 

WflXX  wftXX 

current  1^  (t)  is  computed  by  integrating  (♦, t)  over  the  poloidal 

angle 
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VIII.  Tvo  Applications  for  the  Shell  Model 

When  (b-a)/a  tends  to  zero  but  «(b-a)/a  reaains  finite,  i.e.,  for  the 
shell  Bodel,  Eqs.  (57a),  (57b)  and  (58)  siaplify  considerably,  i.e., 

c  c 

+  tn(l  ♦  -  2  cos{^-a(t))]  (86a) 

A  ft 

8r  * 

*  -r  -  2) 

inside  the  ring  (p'  <  r^), 

aJS  (p,f,t)  -  ^  I^(t)  [2tn  |r  (86b) 

.  U  (i  .  -  2  iHp-  cos(6  -  «(t)))] 

ft  ft 

8r  * 

.  (tn  -  2)  U^<t)  *Y  S  ©*  (“«'<'>  **”♦) 

Bal 

outside  the  ring  but  inside  the  torus  (p  <  a),  and 

A^‘(p,*,t)  .  (tn  2)  U„„(t)  (87) 

m 

*  Z  ■  (t)"  '*"■*)’ 

Bal 
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outside  the  torus  (p  ^  a).  Notice  that  although  t)  is  zero  at  t  - 

0,  its  partial  time  derivative  is  not  zero  at  t  ■  0,  because,  as  explained 
in  the  previous  section,  the  fast  diffusing  terms  which  render  the  electric 
field  zero  at  t  -  0  outside  the  torus,  are  missing  in  the  shell  model. 

Also,  notice  that  at  t>o  the  vector  potential  is  the  sum  of  the  external 
vector  potential  of  the  ring  and  its  image,  but  for  only  the  vector 

potential  of  the  ring  remains  present. 

As  a  second  application,  consider  the  case  in  which  the  ring  moves  on 
a  circle,  i.e.,  A  is  constant,  a  -  t»t,  and  its  current  is  a  step  function 
of  time,  i.e.,  I^tt)  -  1^  0(t).  From  Eqs.  (52a)-(52c)  we  have: 


m  cosmut  atCq  sinmttt  -  e 
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(91b) 


,»xm  sinmur  -  wx-,  cosmwt  +  wx-,  e 
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and  the  vector  potential  becomes: 
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-  tn  (l  *  -  2  ^ 

^  sin(^(ot) 

1  -  ^  cos(f-«»t) 
a 

inside  the  ring  (p'  ^  r^), 

{*,*,t)  .  ^Iof2  (tD  2]  (92b) 

-  2^t.n  —g - 2je  +  tn  ~  ^  ^  cos(^-»t)j 

a  a 

-  - - - j  infl  +  -28^  cos(4-ttt)] 


1  +  («tp)' 


Artan 


COS(^ttt) 
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Arctan 


sin(4-ot) 

fl 

1  -  -^^cos(^-«t) 


-  Arctan 
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1  _  e 

.2 


-t/Tp 


outside  the  ring,  but  inside  the  torus  (p  <  a),  and 

A™'(.,*.t)  .  2  (in  2)  (93) 

,  8r  .  /  /a\2  a  \ 
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outside  the  torus  (p  ^  a).  There  are  tvo  extreme  cs^es  o£  interest: 
i)  when  ««Tq  «  1;  then  the  vector  potential  is  the  sane  as  that  of  a 
motionless  ring,  ii)  vhen  mTq  »  1;  then  there  is  diffusion  of  the  terms 
associated  with  the  loop  tine  but  the  image  fields  do  not  dissipate  to 
zero,  but  they  follow  in  phase  the  circular  motion  of  the  ring. 
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IX< - . Conclusion 

The  diffusion  of  the  self  magnetic  field  of  a  beam  inside  a  toroidal 
conductor  is  governed  by  three  different  diffusion  tines.  The  loop  tine 

*^00  is  responsible  for  the  diffusion  of  the  fields  into  the  hole  of  the 
doughnut  and,  after  a  fev  e-folds  of  the  fast  diffusion  time,  the  tine 

behavior  of  the  vail  current  is  completely  determined  by  The 

"cylinder"  diffusion  time  x^  is  responsible  for  the  dissipation  of  the 

image  fields  which  are  present  initially,  but  they  vanish  after  a  fev 

e-folds  of  1^,  if  the  ring  current  does  not  vary  with  tine.  Finally,  the 

fast  diffusion  time  is  responsible  for  the  electric  field  outside  the 

conductor  to  be  zero  initially,  but  it  acquires  approximately  the  value 

associated  with  the  shell  model  after  a  fev  e-folds  of  Xpj^. 

After  a  fev  e- folds  of  the  loop  time  and  if  the  ring  current  does  not 
vary  with  time,  the  vector  potential  becomes  equal  to  that  in  the  absence 
of  the  conductor.  In  addition,  to  zero  order  in  the  toroidal  corrections, 
the  radial  component  of  the  self-magnetic  field,  which  is  responsible  for 
the  beam  trapping,  is  independent  of  the  loop  time.  Therefore,  the  time 
scale  of  the  trapping  mechanism  should  be  independent  of  x^^.  Reliable 
results  close  to  the  conducting  vail  can  be  obtained  only  by  numerical 
computation  of  the  poles  and  by  including  a  very  large  number  of  terms  in 
the  series  expansions  of  the  vector  potential  and  the  fields.  But  in  the 
extreme  case  of  the  shell  model,  the  results  are  exact  everywhere  inside  as 
veil  as  outside  and  in  the  vicinity  of  the  torus.  This  model  provides 
quite  an  accurate  description  of  the  diffusion  process  for  a  toroidal 
conductor  with  a  thin  vail,  except  during  the  first  fev  e-folds  of  Xp^^, 
since  the  effect  of  the  fast  diffusing  terms  is  not  included.  Due  to  the 
simplicity  of  the  shell  model,  it  is  rather  easy  to  compute  the  first  order 
toroidal  corrections.  These  results  will  be  reported  in  a  future 
publication. 
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Appendix 

An  expression  for  f^(a),  given  by  Eq.  (54b),  will  be  derived  here, 
correct  to  second  order  in  e/x^.  Using  the  eultiplication  theoree  for  the 
Bessel  functions,  ve  obtain  the  relation 

I  '.<«>  •  <*» 

O  1 


.)[• 


B-l+k 


<’■«>  -  Vl<-o>Vuic«*o>]- 


With  the  help  of  the  identities 


W  -  ^<-0)  -  -  s;- 
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and  for  I  >  1,  2,  ...  we  can  also  show  that 

.^[.2124^.0^-1,]. 

o 

'^m+l^’^o^  "  \+e^*o^  '^in+4t^*o^ 


(A2a) 

(A2b) 


(A3a) 

(A3b) 

(A3c) 


(A4a) 


55 


_  4t(2t-l)(« 


(A4b) 


^■+4e+l^*o^  ■  ‘^■♦4t+l^*o^ 


2  rAeim+2t±2J 


[■ 

n*-.  L 


b  >•  *o 


W]. 


(A4c) 


‘^■+1^*0^  ^■+4£+2^*o^  ■  '^■+4£+2^*o^ 


2Ul)(m^^2fc^•l)(ln^^2fc+2^ 


<A4d) 


Therefore,  ve  conclude  that 


•al"  1 


Ir  «■<•)  ■  ^  -  »ln  f 


2  2 


-ill 


.k  (k-fl)(B->-k-»2 
'  (2k+4)! 


4  V  t(2e-l)(m»2f)(m-^2t-i^l)  fb 

,  2  (4t+i)i  r 


'o  £-1 


2-a2  1 

17- 


«  9  9 

4  V  fc(4£-t-l)(iii+2£-.-l)(m.^2£-t-2)  fb^  -  a^  ,  ^ 
^  2  A  (4t+3)I  I  2a^  *oJ 

o  £»1 


where  terms  of  order  1/x^  and  higher  have  been  omitted.  A  straightforward 
and  lengthy  calculation  leads  to  the  relations 


Z.  -.k  (k-»l)(m->-k>2 
(2k+4)! 


(A6a) 


56 


■  ■  -  j  -  l^jsln  z  -  jli+l— jcos  z, 
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Substituting  Eqs.  (A6a),  (A6b)  into  Eq.  (A5),  ve  conclude  that 
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Ve  8««  that  the  zaroas  of  f^(a)  are  given  by  the  rt'  on 


2 

correct  to  order  (n/x^)  .  If  (b  >  a)/a  «1,  and  x^  »1,  the  relation  above 
siaplifies  to 
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Table  I.  Parameters  o£  the  run 

Torus  major  radius  r 

o 

Torus  inner  minor  radius  a 
Torus  outer  minor  radius  b 
Strong  focusing  radius 
Strong  focusing  current 
Strong  focusing  Periodicity  m 
Vertical  magnetic  field 

zo 

Toroidal  magnetic  field  ®eo 
Beam  relativistic  factor  r 
Beam  minor  radius  r 

c 

Beam  current  I. 

c 

Vail  resistivity 
Intermediate  frequency 


shown  in  Fig.  2 

100  cm 
15.2  cm 
15.217  cm 
23.4  cm 
24  kA 
6 

26  Gauss 
4000  Gauss 
1.69714 

0.5  kA 
8.84  mS  cm 

9 

1.8  X  10  sec 
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2 


Electron  Beam 


Fig.  1  System  of  coordinates 
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Time  (/isec) 


Fig.  2  Beam  centroid  orbit  [(a)  and  (b)]  and  relativistic  factor  vs.  time 
[(c)]  from  the  numerical  integration  of  the  ring  equations  of  motion 
coupled  with  the  diffusion  fields. 


